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ABSTRACT 

Cosmological black holes (CBH), i.e. black holes with masses larger than 10 1 4 Mq, 
have been proposed as possible progenitors of galaxy voids ([Stornaiolo 20021) . The 
presence of a CBH in the central regions of a void should induce significant gravita- 
tional lensing effects and in this paper we discuss such gravitational signatures using 
simulated data. These signatures may be summarized as follows: i) a blind spot in 
the projected position of the CBH where no objects can be detected; ii) an excess of 
faint secondary images; iii) an excess of double images having a characteristic angular 
separation. All these signatures are shown to be detectable in future deep surveys. 
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1 INTRODUCTION 

Voids are among the largest structures known in the Uni- 
verse with typical diameters ranging between 20 and 85 
Mpc. Among the various models so far proposed for their 
form ation, we may quote th e one originally introduced 
in (|Friedmann and Piran 200 lh . accordingly to which voids 
form from the evolution of negative primordial perturba- 
tions in the density field. More in detail, in this model void 
formation is the result of two correlated processes. The first 
one is the comoving expansion of these negative fluctua- 
tions. The second arises from the biased galaxy formation 
picture: galaxies are less likely to form in the underdense 
regions created by this expansion. Several N-body simula- 
tions of this formation mechanism base d on the cold dark 
matter scenario (cf. ((Benson et al 20 03)) produced results 
which are consistent with observational data. However, it 
needs to be stressed that the current observational samples 
of data ( void surveys, mor phological classification of void 
galaxies) (Ro ias et al 2 005) are too small to allow any con- 
clusive test of this or other models. 

Anot her formation mechanism was proposed by Stor- 
naiolo (IStornaiolo 2002T) . In this scenario the collapse of ex- 
tremely large wavelength positive perturbations led to the 
formation of low density/high mass black holes (Cosmolog- 
ical Black Holes or CBH) . Voids are then formed by the co- 
moving expansion of the matter surrounding the collapsed 
perturbation. This model implies that at the center of voids 
should still exist a very massive (M > 1O 14 M ) CBH, 
which should be detectable through the gravitational lens- 
ing effects induced on background galaxies. In this paper we 



discuss what these effects are and whether they might be 
actually observable or not. 

This paper is structured as follows. In Section we 
present the simulations which were performed in order to 
derive the possible gravitational signatures of the CBH, sum- 
marized in Sectional In Section 0] we discuss whether such 
effects may or may not be observed in existing or ongoing 
surveys. In Section we draw some conclusions. In a subse- 
quent paper we shall discuss the weak lensing effects possibly 
induced by a CBH on a background galaxy distribution. 



2 THE SIMULATIONS 

In order to evaluate the gravitational lensing effects induced 
by a CBH located in the center of a void we produced two 
sets of simulations: one set for a reference unperturbed uni- 
verse, i.e. for a void in an otherwise uniform universe (with- 
out the CBH) and a second set obtained from the previous 
by adding a CBH in the center of the void. 

In both cases we assumed, for simplicity, an Einstein-de 
Sitter cosmological background model, with Qm = 1, Ho = 
lOO/i Km s _1 Mpc -1 and no cosmological constant. We wish 
to stress that the only effect of a non zero cosmological con- 
stant would be to change the relation D = D(z,Ho, Qm, Ha) 
between the distance D and the redshift z of the background 
galaxies which is involved in the calculation of the apparent 
magnitudes and of the angular positions. U sing approximate 
expressions ([Kantowski and Thomas 200 lh it is easy to see 
that over a redshift range [0,0.4] the maximum deviation 



2 M. Serpico 



■<*5r . 

C 2 

m 



o 



B' 



js • mi 

C CBH 



A" 



Figure 1. Schematic layout of how the CBH gravitational effects 
affect the formation of double images as a function of the angular 
distance of the source from the CBH (see text for a more detailed 
explanation). 



from the standard Einstein-de Sitter model would be of the 
order of 10%. 

As mentioned, the reference universe is described by an 
Einstein-de Sitter background containing a spherical void of 
radius R VO id located at the comoving distance D voi d from the 
observer which does not produces any detectable deflection 
or magni fication effects on b ackground galaxies (as pointed 
out in (|Amendola et al 19991) ). 

The perturbed universe is instead described by a Swiss- 
Cheese model, with a CBH in the center of the void, which, 
as above, has comoving radius R VO id and is located at the 
comoving distance D VO id from the observer. The CBH has a 
mass 



4 3 

M = - TV Q cb h pcritRyoid 



(i) 



where Q c bh is the density parameter of the CBH. It 
needs to be stressed that while the background universe 
is described by a FLRW metric, inside the void region 
holds a Schwarzschild metric. The boundary conditions 
at the transition between the two regimes are discussed 
m (iKantowski 1969ft . We also assume that changes in the 
sources redshift, due to the propagation of light across the 
edge of the void, is negligible. In other words, we describe 
the void-CBH system as a Schwarzschild lens enclosed in a 
Einstein-de Sitter cosmological model. 

More in particular we simulated a slice of universe covering a 
solid angle defined by the void diameter and a redshift range 
comprised between the far edge of the void and z = 0.4. This 
conical volume was populated with a randomly distributed 
(in the comoving frame) galaxy population drawn from the 
r-band SDSS (Slo an Digital Sky Sur vey) luminosity func- 
tion for the field (Bla nton et al 2003F ). Evolutionary effects 
induced by the redshift on the luminosity function were ne- 
glected. Furthermore we treated galaxies as material points 
since their physical extension is not relevant to the following 
discussion. 

induced by the CBH were 
weak lensing approximation 
assuming that the light from 
distance from the CBH much 
larger than its Schwarzschild radius (about 30 arcsec for 
M = 10 14 M© and D void = 50 Mpc). 



The lensing effects 
then derived in the 
(ISchneider et al 19921) by 
the sources passed at a 



The deflection angle produced by the CBH for radiation 
approaching with impact parameter £ is given by 

A 4GM 



c 2 £ 



(2) 



For a source at distance D s from the observer, the Ein- 
stein angle can be written in the form: 



a 



4GM D ds 

D vo idD S 



(3) 



where Dds is the distance of the source from the lens. 

If j3 is the unlensed angular position of the source with 
respect to the observer, we know that the effect of the lens 
will be the creation of two images of the source with angular 
positions: 



and with magnifications given by: 



Ml,2 = i 



v//3 2 +4 



+ V£±i ± 2 



(4) 



(5) 



where (3 = (3/ao. 

Simulations were then performed for different values of 
Qcbh and R V oid (0.05 to 0.4 with step 0.05 and from 10 to 
20 Mpc, step 2, respectively). 



3 QUALITATIVE DESCRIPTION OF 
OBSERVABLE QUANTITIES 

The simulations showed that the CBH leaves three different 
types of signatures on the background galaxy distribution. 
In order to better quantify what happens we refer to Fig. 1 
which shows the images produced in three different relative 
positions of source and lens. When the "real" object (A) is 
very close to the CBH, we have the formation of two images, 
namely A' and A" which are respectively outside and inside 
the Einstein Radius and very close to it. Both images are 
brighter than the unlensed image and A' is brighter than A" . 
Then when the source moves away from the CBH, the am- 
plification factor with respect to the secondary image tends 
to 1. 

The dashed inner circle marks the position of such locus. 
If the source B lays on this circle, then it will produce two 
images B' and B" . B' is outside of the Einstein Radius at 
a larger distance than A' and is brighter than B, while B" 
lays inside the Einstein Radius closer to the CBH than A" . 
Finally let us consider the case of a source C which, if un- 
perturbed would fall outside of the Einstein Radius. Also in 
this case we shall see two images C' and C" . The first one 
will almost coincide with the position of C and will have an 
almost identical brightness, while C" will be almost invisible 
and very close to the CBH. 

When a random distribution of background galaxies is 
considered, the overall effect of the lensing will be the for- 
mation of 4 different areas on the sky, namely the regions 
A, B, C and D, shown in Fig. 2. 

The inner circle A, which we call 'blind spot', is the region 
where no image can be detected due to the increasing de- 
magnification of secondary images when they move towards 
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^CBH /Rvoid 


10 


12 


20 


0.05 


16±5 


25±7 


36±8 


0.2 


33±8 


36±9 


53±10 


0.4 


39±8 


44±9 


69±12 



Table 1. Average value of Obiind (in arcsec) as a function of 
&CBH an d R V oid (i n Mpc). The quoted errors are the r.m.s. of 
the 1000 individual simulations. 

the center of the void. The size of the blind spot depends on 
the mass of the CBH and on the density of the background 
galaxy distribution. 

The annulus B is characterized by the presence of a large 
number of secondary images, and it is the second observable 
feature associated with the CBH. 

The third zone C, is an annular zone that we call the 'deficit 
zone'. It encompasses the average Einstein Radius of the 
galaxy sample and is characterized by a relatively low num- 
ber of background galaxy images. It can be understood re- 
minding that, at this angular distance from the center of the 
void, one can find only those primary and secondary images 
which originates from sources having angular position falling 
well within the Einstein Radius. Therefore, in this annulus 
it is expected to find a lack of both primary and secondary 
images: as the source moves away from the Einstein Radius, 
the primary image tends to coincide with the original posi- 
tion of the source while the secondary image becomes fainter 
and moves more and more towards the CBH. 
Finally, the D zone does not present any particular feature 
produced by lensing, and coincides with the homogeneous 
background galaxy distribution. 

In total we run 48 sets of simulations assuming the void 
distance at 50 Mpc (i.e. matching the distance of the nearest 
void) and covering a grid denned by: Qcbh = 0.05 — > 0.4 
with step 0.05, and R VO id = 10 — > 20 Mpc with step 2 
Mpc. For each grid point the procedure was iterated 1000 
times randomly changing at each iteration the positions of 
the background galaxies. Each simulation produced a cata- 
logue of galaxy positions and magnitudes and each group of 
simulations was then used to derive average quantities. 



4 RESULTS 

The main observational signatures left by the CBH can be 
summarized as follows. 




Figure 2. The three regions described in the text. Region A 
(greatly enlarged to make it visible) is the blind spot; Region B: 
region where we expect the excess of secondary images; Region 
C: deficit region (the circle inside the C region shows the average 
Einstein Radius). 
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Figure 3. The figure shows in red the average number counts as 
a function of the apparent magnitude in the r band, obtained, 
respectively, in the outer region (see text) and in black those 
obtained in the inner region. 

size is of the same order of the average angular separation 
between bright galaxies at intermediate redshift. Possible 
effects connected with the distortion induced on extended 
background objects which happen to fall near the line of 
sight of the CBH will be addressed in a forthcoming paper. 



4.1 Blind spot 

For every simulation, we determined an estimate of the an- 
gular radius Ouind of the blind spot, defined as the minimum 
angular distance from the CBH of the secondary images 
(brighter than mu m = 23.5) and then, at each simulation 
grid point we took the average value over the 1000 simula- 
tions. In Table 1 we list some representative values. 

As expected, the size of the blind spot, increases with 
Qcbh and with R VO id- Assuming, for instance, a typical 
value of 30 arcsec (cf. Table 1), it is apparent that the blind 
spot should be rather difficult to observe. In fact, while it 
should be easily detectable in very deep number counts, its 



4.2 Galaxy number counts and radial profile 

The presence of the CBH affects the number counts. In or- 
der to estimate the size of such effect and in absence of a 
priori information on the size of the void we adopted the 
following procedure. First we introduced an annular zone 
(defining an inner and outer region) centered on the blind 
spot. The radius of the zone was then found by maximiz- 
ing the difference between the average galaxy counts in the 
inner and outer regions. 

As it can be seen in Fig. 3, the average number counts 
associated with the inner area show a systematic differ- 
ence with respect to those extracted from the background. 
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This effect however becomes significant only at magnitudes 
fainter than ~ 21.0. 

In order to quantify such difference we performed a 
Kolmogorov-Smirnov test on both distributions of points. 
Each data set was split into two parts including galaxies 
brighter or fainter than 21 mag, respectively. The brighter 
parts of the distribution do not present any statistically sig- 
nificant difference, while for the fainter parts we derived a 
probability higher than 98%, that the two samples are drawn 
from different populations. 

As it was discussed in the previous paragraphs, the 
presence of a CBH induces a typical pattern in the num- 
ber counts radial profiles. Such pattern is characterized by 
a peak in the range of distances intermediate between the 
blind spot radius and the average Einstein Radius (caused 
by the secondary images concentration), followed by a dip, 
which corresponds to a slight underdensity of objects and 
than at distances comparable with the Einstein Radius it 
raises up again to smoothly reach the value expected for the 
background galaxy distribution. 

In Fig. 4a we show the number counts profile extracted 
from the simulation grid-point at Qcbh = 0.2 and R VO id = 
12 Mpc. The first point of the profile is located in the blind 
spot and is followed by an isolated peak which rapidly falls 
in the dip. 

The only two other examples of possible radial profiles 
obtained when not placed on the center of the void are il- 
lustrated in Figg. 4b and 4c. The first image, which can be 
roughly described by an initial peak, larger than the previ- 
ous case, located at the first steps of the profiles followed by 
the dip, is the typical pattern generated when the central 
point of the radial profile is positioned inside the secondary 
images overdensity; the second profile, characterized by the 
presence of an initial low spot followed by two distinct peaks 
and the dip, is reproduced when the center of the radial den- 
sity profile lays out of the circular overdensity created by the 
secondary images, and inside the deficit annulus. 



4.3 Multiple images 

The third signature comes just from the pairs of double im- 
ages produced by the CBH. We expect angular separations 
for these pairs of the order of two times the average Ein- 
stein's angle of the sample. This means angular separations 
of the order of some arcminutes. This forecast is particularly 
striking because if these doubles really exist, the only way 
to observe them practically is to have an estimate of their 
wide angular separation. We measured the distribution of 
the angular separation between the double images produced 
by the CBH obtaining the two points angular correlation 
function w(0) for th e simulated galaxy distri bution. A stan- 
dard Landy-Szalay (|Landv and Szalav 1993i) estimator was 
used: 

(DD) + (RR) - 2(DR) 



w{0) 



(RR) 



(6) 



where (DD), (DR), (RR), are pair counts in bins of 6 ± SO 
of: data-data, data-random and random-random points, re- 
spectively. The statistic has been demonstrated to be close 
to a minimum of variance estimator an d to be robust with re - 
spect to the number of random points (|Kerscher et al 2 000). 
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Figure 4. Galaxy number counts radial profiles obtained by inte- 
grating over annular regions centered in different regions. Panel a 
- centered on the blind spot (zone A); Panel b - center falls in the 
secondary images overdensity region (zone B); Panel c - center 
falls in the deficit region (zone C). 



In Fig. 5 we show the above defined correlation function ob- 
tained for the simulation grid-point at Qcbh — 0.2 and 
Rvoid — 12 Mpc; a well defined peak corresponding to an 
angular distance of 6 arcmin, comparable with the average 
diameter of the secondary images overdensity region (Zone 
B), is visible. Moreover, a slight anticorrelation is found at 
distance greater than 10 arcmin, while for a random dis- 
tribution of points no correlation of any sort should be de- 
tected. It needs to be stressed that the peak observed at very 
short angular separations is an artifact produced by the fact 
that in our simulations galaxies are assumed to be point-like 
and would disappear if galaxies were approximated with ex- 
tended objects. 
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Figure 5. The angular 2-point correlation function derived from 
the simulated galaxy distribution (see text). 

5 DISCUSSION AND CONCLUSION 

The results of our simulations which were based on a min- 
imal set of assumptions and can therefore be regarded as 
quite general and robust in their final predictions, clearly 
show that the presence of a CBH close to the center of a 
void would leave unmistakable signatures on the background 
galaxy distribution as a result of the gravitational lensing 
properties of the CBH. Unfortunately, they also show that 
such signatures can be detected only at faint light levels, i.e. 
at magnitudes fainter than the completeness limit of most 
existing photometric surveys which include voids in their 
field. Furthermore, the only survey which, at least in the- 
ory, should be deep enough to allow at least the radial profile 
test above described (nam ely the Sloan Digital Sky Survey, 
cf. (jStoughton et al 20021 )1 does not satisfactorily cover any 
previously known void. It needs to be stressed however that 
our results clearly show that such tests will be possible on 
any of the planned deep digital surveys which will become 
available in the near future (cf. for instance the VST extra- 
galactic survey). 
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